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C ■ Abstract. We present our theoretical results on an effective two-band double-exchange model 

^ on a pyrochlore lattice for understanding intricate phase competition in Mo pyrochlore oxides. 

The model includes the twofold degeneracy of e' g orbitals under trigonal field splitting, the 
interorbital Coulomb repulsion, the Hund's-rule coupling between itinerant e' g electrons and 
localized oi B spins, and the superexchange antiferromagnetic interaction between the a\ g spins. 
By Monte Carlo simulation with treating the Coulomb repulsion at an unrestricted-type mean- 
field level, we obtain the low-temperature phase diagram as functions of the Coulomb repulsion 
R , and the superexchange interaction. The results include four dominant phases with characteristic 

spin and orbital orders and the metal-insulator transitions among them. The insulating region 
is characterized by a 'ferro'-type orbital ordering of the e' g orbitals along the local (fit) axis, 
irrespective of the spin ordering. 
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1. Introduction 

Mo pyrochlore oxides, whose chemical formula is generally given by -R2M02O7 (R is rare earth), 
are known to exhibit interesting phase competition and related phenomena. By chemical 
substitution of R cations, the system changes from a ferromagnetic metal (FM) to a spin-glass 
(SG) insulator systematically depending on the ionic radius of R cation pQ. On the other hand, 
in applied external pressure, FM collapses and turns into a peculiar paramagnetic metal (PM) 
with showing an intermediate SG metallic state, whereas the SG insulator remains rather robust 
against pressure [2J. This variety of phases and competition among them offer an opportunity to 
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systematically study many fascinating phenomena in strongly correlated electron systems, such 
as metal-insulator transition and anomalous transport properties [3]. 

The basic lattice structure of -R2M02O7 is composed of two intervening pyrochlore lattices: 
One is formed by Mo cations and the other is by R. The pyrochlore lattice is a three-dimensional 
network of corner-sharing tetrahedra, which is known to possess strong geometrical frustration. 
Mo cation is nominally tetravalent and has two electrons on average. The crystal field from 
trigonal distortion of M0O6 octahedra splits three-fold ti a levels of Mo Ad electrons into a\ g 
singlet and e' doublet. The trigonal axis is along the direction toward the center of each Mo 
tetrahedron, which is called the local (111) axis. With aligning their spins parallel according 
to the Hund's rule, one of two electrons occupies the lower a\ g level, and the other comes in 
the doubly-degenerate e' g levels. (See, for example, Fig. 1 in Ref. [2].) Hence spin, charge, and 



orbital degrees of freedom are all active in these compounds, and the frustration of the lattice 
structure promotes the competition among them by preventing a simple ordering. In addition, 
the coupling between Mo Ad electrons and -R-site rare-earth moments brings about further 
complications. In the following, however, we will focus on the phenomena commonly observed 
regardless of magnetic or nonmagnetic R cations, such as the phase competition mentioned 
above, and neglect the R sites: We will concentrate on the Mo pyrochlore network alone and 
the intrinsic physics from Mo Ad electrons in this paper. 

Several theoretical studies have been done for understanding the phase competition in the 
Mo pyrochlore oxides. A pioneering work was done for the chemical substitution at ambient 
pressure by Solovyev by using the LDA+C7 calculation [3]. The calculated electronic structure 
indicates that the lower a\ g level has rather localized nature with a narrow band, while the e' 
levels form a relatively wide band and bears itinerant nature. The observation suggests that 
the system can be viewed as a double-exchange (DE) system, in which localized a\ g spins are 
coupled with itinerant e' electrons via the ferromagnetic Hund's-rule coupling. This explains 
well the emergence of FM in the compounds with large ionic radius of R, such as i?=Nd and 
Sm. The calculation while changing the value of U shows that the localization from FM to SG 
insulator by the substitution is basically understood by the correlation effect (or the bandwidth 
control). These calculations, however, were performed with assuming a simple collinear ferro or 
antiferromagnetic (AF) ordering; hence, the detailed nature of the frustrated insulating phase, 
such as the SG behavior and the orbital state, was not clarified. 

Recently, the authors have studied the phase competition under applied pressure El [7] . 
The results show that the change from FM, SG metal, to PM is well reproduced by Monte Carlo 
(MC) calculations for an extended DE model including the superexchange (SE) AF coupling 
between localized spins. In particular, the peculiar diffusive nature of PM is explained by the 
strongly incoherent state originating from competition between the DE ferromagnetic interaction 
and the SE AF interaction [5]. The SG metallic state is also accounted by an instability toward 
phase separation related to the competition [7J. The agreement implies that the physical pressure 
causes a different effect from the chemical pressure (chemical substitution), as implied in the 
experimental report [2]: The former tends to increase the SE AF coupling between the a\ g 
spins dominantly, while the latter mainly modifies the bandwidth of itinerant e' electrons. MC 
calculations so far were performed for a single-band model with neglecting the Coulomb repulsion 
between itinerant electrons; for comprehensive understanding of the phase competition in the Mo 
pyrochlore oxides, it is necessary to extend the analysis to incorporate the electron correlation 
as well as the e' g orbital degree of freedom. 

In this contribution, we investigate the effect of electron correlation and orbital degeneracy 
on the phase competition in the extended DE model. The model explicitly includes the twofold 
orbital degeneracy and the Coulomb repulsion between itinerant electrons. Thermodynamic 
properties of the model are studied by MC calculations, in which the electron correlation is 
handled by an unrestricted-type Hartree approximation. The results indicate that, in addition 
to the phase transition from FM to PM previously obtained in the single-band model, a metal- 
insulator transition is caused by the electronic correlation accompanied by a 'ferro'-type orbital 
ordering. We discuss the phase diagram in relation to the distinct behaviors in experiments for 
the chemical substitution and the applied pressure in the Mo pyrochlore oxides. 

2. Model and method 

As an extension of the previous studies El [7] , we here consider the DE model with twofold 
orbital degeneracy. The Hamiltonian is given by 

7* = EEEW( C L C J^ + H - C ')+ 2^ 5Z Y, U <xP<x'P' C ia>r C hcT> c iP'cr'Ci*><T 
(ij) a/3 o i afia'fi' aa 1 
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where c iacr (ci aa ) creates (annihilates) an electron with orbital a and spin a at site i, a denotes the 
Pauli matrix, and Si represents the localized spin at site i. The first term describes the electron 
hopping between nearest-neighbor sites (ij) on the pyrochlore lattice; the transfer integrals 
depend on the orbital indices a and (5 which represent two e' orbitals. The second term denotes 
the onsite Coulomb interactions between the e' electrons, including the intra and interorbital 
repulsions, the Hund's-rule coupling, and the pair hopping. The third term denotes the Hund's- 
rule coupling between itinerant electrons and localized spins Si. The last term is the SE AF 
coupling between the neighboring localized spins. A similar model has been studied on the 
square or cubic lattices in the context of colossal magnetiresistive (CMR) manganites [8]. 

Our purpose is to clarify thermodynamics of the model given by Eq. ([1]) and to elucidate 
the competing phases including correlated disordered states, such as the diffusive PM or the 
Mott insulating state. For this purpose, some sophisticated method beyond a simple mean- 
field approximation is necessary. The model ([T]), however, is very complicated and difficult to 
handle beyond a simple mean-field approximation. Hence, we introduce several simplifications 
as follows. First we take the limit of large Jr. In this limit, as often considered in the study 
of CMR, spins of itinerant electrons are aligned parallel to the localized spin at each site, and 
the transfer integrals are modulated by the relative angle between the neighboring localized 
spins. At the same time, among may contributions in the second term in Eq. (TjQ), only the 
interorbital interaction between the same spin electrons remains relevant. Furthermore, we 
consider the localized spins Si as classical vectors with the renormalized length \Si\ = 1. The 
simplifications were shown to give reasonable results in the study of a single-band model for the 
phase competition under pressure in the wide-bandwidth compounds [3 El [7]. Consequently, by 
taking the local spin axis along the direction of localized spin at each site, the model can be 
described by spinless fermions in the form 

H = X1X1 iiaJpiCiaCjP + H - c + "9 5UI «ia"i/3 + ^AF Yl ^ ' &3> ( 2 ) 
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where di a is the spinless fermion operator and hi a = cl a 5i a . Here, U is considered as an effective 
interorbital interaction, which controls the double occupancy in the e' orbital with the same spin 
electrons under the strong Hund's-rule coupling; its origin is traced back to the e' interorbital 
Coulomb repulsion and the Hund's-rule couplings between e' g -e' g and e' g -ai g orbitals in Eq. ([1]). 
Hereafter, we consider the filling with one electron per site on average. 
In Eq. ([2]), the effective transfer integral is given by 

iiaj/3 = t a i3 cos -£ cos -^ + sin -£ sin -|- exp{-i(<fo - <f>j)} , (3) 

where (Oi,<pi) denotes the angle of Si. To estimate the value of t a p, it is necessary to take 
account of the anisotropy of the e' orbitals and the relative angle of Mo-O-Mo bond. Following 
the calculation in Ref. [9], we estimate the orbital diagonal and off-diagonal elements: it turns 
out that the ratio between them depends solely on the Mo-O-Mo bond angle 5 as 

tg^/3 _ 3 - cos 5 
t a= /3 3 + cos 5 

Thus, the orbital off-diagonal element is larger than the diagonal one for 5 > 90°. In Mo 
pyrochlore oxides, 5 is estimated around 130° at ambient pressure [10]. In the following 



calculations, for simplicity, we consider only the dominant off-diagonal element and fix t a ^p = — 1 
and t a= p = 0, although the transfer integrals might change depending on both chemical and 
physical pressure. 

We study the model ([2]) by MC simulation in a similar way as in the previous studies 
for the single-band model [5j El \7\. To facilitate the MC calculations without negative sign 
problem, we introduce an unresticted-type mean-field decoupling to the Coulomb interaction 
in the second term in Eq. ([2]). Here, as the simplest approximation, we employ the Hartree 
approximation; hi a riip — > (ni a )nip + ni a {fiip) — (nj Q )(nj / g). The MC sampling is performed to 
update the configurations of localized spins {Si}. In each MC step, the mean fields {(rii a )} are 
recalculated by using the eigenvalues and eigenvectors obtained from the exact diagonalization 
of the Hamiltonian (|2|) . We confirmed the convergence of the mean fields within the statistical 
errorbars after thermalization. We also employ the averaging over the twisted boundary 
conditions to suppress the finite size effect [6j. In the following, we show the result for 16 
site system with averaging over 8 3 different twisted boundary conditions. Typically, we perform 
2000 MC steps for measurement after 1000 MC steps for thermalization. 



3. Result 

Figure [1] shows the MC results for the model ([2} at a low temperature T = 0.04 while changing 
the Coulomb repulsion U and the SE AF coupling Jaf- We compute the square of magnetization 
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the charge compressibility Xc and the z component of orbital polarization t z , by 
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respectively. Here the bracket (• • •} denotes the thermal average, N s is the number of sites, \i is 
the chemical potential, and n is the electron density. 
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Figure 1. MC results for the model given by Eq. (J2} at (h) = 1 and T = 0.04: (a) square of 
the magnetization m? [Eq. (JSJ)], (b) charge compressibility Xc [Eq. ©], and (c) z component of 
orbital polarization t z [Eq. ((7J)]. The results are obtained for 16 site cluster by taking average 
over twisted boundary conditions with 8 3 different fluxes. 



As shown in Fig. [TJa), in the small Jaf region, the magnetization is large, indicating that the 
system is ferromagnetic. The ferromagnetism collapses as Jaf increases. On the other hand, 



the magnetism is relatively insensitive to U. For small U, the charge compressibility Xc is finite 
as shown in Fig. QJb), indicating that the system is metallic. Hence, the system is FM for small 
Jaf, while it turns into PM as Jaf increases in the small U region. This behavior is qualitatively 
the same as observed in the single-band model |5j. In the single-band case, phase separation 
appears between the two phases [7]; we expect a similar behavior in the present two-band case, 
but it is not clearly seen in the present calculation presumably because of the small size cluster. 

As U increases, Xc decreases rapidly to zero, as shown in Fig. \H[b). This signals the 
metal-to-insulator transition. The insulating state is characterized by an orbital ordering. 
Figure [T^c) shows that the orbital polarization becomes nonzero almost simultaneously with the 
disappearance of the charge compressibility. This means that the insulating state shows a 'ferro'- 
type orbital order by choosing one of two e' orbitals. The reason for the 'ferro'-type orbital order 
is understood from the strong coupling picture as follows. In the large U limit, the second order 
perturbation in the transfer integrals leads to the so-called Kugel-Khomskii-type spin-orbital 
exchange couplings [IT]. In the present model, as seen in Fig. [JJa), the spin configuration is 
governed by the competition between DE and SE interactions, not by U. Then, the 'ferro'- 
type orbital configuration is favored via the spin-orbital exchange couplings, irrespective of spin 
configurations, because the dominant hopping processes in the perturbation are the orbital off- 
diagonal ones by t a ^p in the present model. The situation is contrastive to the 'myth' for 
the Kugel-Khomskii-type interactions, in which the system usually favors complementary spin- 
orbital configurations, such as spin-AF orbital-ferro or spin-ferro orbital- AF. 

It is worthy to note that the 'ferro'-type orbital order is not a uniform orbital order. The 
trigonal distortion is along the local (111) axis at each site, and hence, the e' orbital polarization 
takes place along the site-dependent local axis. In other words, the 'ferro'-type orbital order is a 
q = four-sublattice order. Thus it does not break the cubic symmetry of the pyrochlore lattice 
structure. 

Summarizing the results, we obtain four dominant phases at low temperatures in the plane 
of U and Jaf ; as schematically shown in Fig. [2j The spin-ferro orbital-para metal in the small 
U and small Jaf region is stabilized by the DE mechanism. The spin-para orbital-para metal in 
the small U and large Jaf region is induced by the competition between the DE ferromagnetic 
interaction and the SE AF coupling under the geometrical frustration. In these two phases, the 
orbital degree of freedom is less important, and hence, the physics is common to the single-band 
case El [7]. On the other hand, the spin-ferro and orbital-'ferro' insulator in the large U and 
small Jaf region is a correlation-driven insulator with development of the orbital ordering. The 
spin-para orbital-'ferro' insulator is driven by the increase of Jaf and the spin frustration. 



spin-para 

orbital-para 

metal 
k 



spin-ferro 

orbital-para 

metal 



spin-para 

orbital-'ferro' 

insulator 



spin-ferro 

orbital-'ferro' 

insulator 



U 



Figure 2. Schematic phase diagram at low temperatures 
suggested from the results in Fig. [1] The dotted (dot- 
dashed) arrow in the figure indicates a possible parameter 
change while changing the physical (chemical) pressure in 
the Mo pyrochlore oxides. See the text for details. 



4. Discussion and concluding remark 

In Mo pyrochlore oxides, as mentioned in Sec. 1, the chemical substitution leads to the transition 
from FM to SG insulator, while the pressure induces the transition from FM to PM. The latter 



is qualitatively explained by the change of Jaf in the single-band model [5], and the present 
calculation for the extended two-band model also gives similar results in the small U regime; 
the parameter change is, e.g., schematically shown by the dotted arrow in Fig. [2j On the 
other hand, the metal-insulator transition by the chemical substitution cannot be understood 
by the single-band physics. From our results, we deduce that the transition from FM to SG 
insulator is described by the transition from the spin-ferro orbital-para metal to the spin-para 
orbital-'ferro' insulator, as indicated by the dot-dashed arrow in Fig. [21 In the spin-para orbital- 
'ferro' insulating state, the spin degree of freedom remains paramagnetic because of the strong 
frustration, but short-range AF correlations are expected by the dominant Jaf- The remaining 
frustration might be released by remnant perturbations neglected in the present model. In fact, 
one of the author and his collaborators recently demonstrated that the frustrated spins have 
instability toward a SG state through the coupling to local lattice distortions in the presence of 
bond disorder, reasonably in accord with the experiments in Mo pyrochlore oxides [12]. Thus, our 
results provide a hint for understanding the nature of the frustrated Mott insulating state in the 
Mo pyrochlore oxides. From these observations, we consider that, despite several simplifications 
introduced, our effective two-band model given by Eq. ([2]) can offer a simple starting point 
for comprehensively understanding the intricate physics -R2M02O7. Further analysis in larger 
system sizes is desired to elucidate the thermodynamics as well as transport properties. 

Why and how the two effects, the chemical substitution and the pressure, lead to such different 
behaviors remain unclear. In the present scenario, these effects are attributed to the different 
control of the Coulomb repulsion (or the bandwidth) and the SE AF coupling between a\ g 
spins. The bandwidth is basically determined by the transfer integrals between the itinerant e' 
orbitals, i.e., t a p in our model, on the other hand, the SE AF coupling Jaf is approximately 
set by t^ /U aig , where t ai and U ai are the transfer integral between the a± g orbitals and the 
intraorbital Coulomb repulsion in the a± g orbital, respectively. Both t a p and t ai are modified 
by the chemical or physical pressure via the change of the angle and/or bond length of Mo-O-Mo 
bonds. In order to clarify how they are modified and to distinguish the two effects, it is highly 
desired to study the microscopic change of the lattice parameters by, e.g., the detailed x-ray 
scattering experiment or the first-principles calculations with structural optimization. 
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